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Abstract. Full counting statistics is a fundamentally new concept in quantum trans- 
port. After a review of basic statistics theory, we introduce the powerful Green's func- 
tion approach to full counting statistics. To illustrate the concept we consider a number 
of examples. For generic two-terminal contacts we show how counting statistics eluci- 
dates the common (and different) features of transport between normal and supercon- 
ducting contacts. Finally, we demonstrate how correlations in multi-terminal structures 
are naturally included in the formalism. 

1 Introduction 

The probabilistic interpretation is a fundamental ingredient of quantum me- 
chanics. While the wave function determines the full quantum state a system 
and its evolution in time, observable quantities are related to hermitian oper- 
ators. Expectation values of these operators determine the average value of a 
large number of identical measurements. However, an individual measurement 
yields in general a different result. Applying this idea to a current measurement 
in a quantum conductor, leads directly to the concept of full counting statis- 
tics (FCS): during a given time interval a certain number of charges will pass 
the conductor. To predict the statistical properties of the number of transfered 
charges we need a probability distribution. The theoretical goal is to find this 
distribution. 

Overview fn this article we give an introduction to the field of full counting 
statistics in mesoscopic electron transport. We will concentrate on the powerful 
technique - using Keldysh-Green's functions - which at the same time also is 
based on microscopic theory. To accomplish this goal we will first review concepts 
of basic statistics, which are relevant for counting statistics. In the next section 
we address the microscopic derivation of FCS using Keldysh-Green's functions. 
In the rest of the article we demonstrate the use of counting statistics in a number 
of examples, like two-terminal contacts with normal and superconducting leads, 
diffusive metals and, finally, multi-terminal structures. But first we review briefly 
the development of the field. 

History Full counting statistics has its roots in quantum optics pQ, where the 
number statistics of photons is used, e. g., to characterize coherence proper- 
ties of photon sources. The major step to adopt the concept to mesoscopic 
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electron transport has been undertaken by Levitov and Lesovik Since then 
the theory of FCS of charge transport in mesoscopic conductors has advanced 
substantially, see Refs. [314) . In Ref. it was shown that scattering between 
uncorrelated Fermi leads with probability T is described by a binomial statis- 
tics P{N) = §T N {1 - T) M ~ N . Here, P{N) is the probability, that out of 
M — 2etoV/h independent attempts N charges are transfered. Furthermore, 
Levitov and coworkers studied the counting statistics of diffusive conductors [5], 
time-dependent problems [H] and of a tunnel junction [7J. A theory of full count- 
ing statistics based on the powerful Keldysh-Green's function method was initi- 
ated by Nazarov |Hj . This formulation allows a straightforward generalization to 
systems containing superconductors [9110) and multi-terminal structures [11112) . 
Classical approaches to FCS were recently put forward for Coulomb blockade 
systems [13114) . and, for chaotic cavities based on a stochastic path- integral ap- 
proach PHI- The field of counting statistics in the quantum regime is closely 
related to the fundamental measuring problem of quantum mechanics, which 
has been addressed in a number of works [6I16I17I18I19) . Expressing the FCS 
of charge transport by the counting statistics of photons emitted from the con- 
ductor provides an interesting alternative to classical counting of electrons |20| . 
Counting statistics has been addressed by now for many different phenomena 

• Andreev contacts [21] 

• generic quantum conductors [EESESEH 

• adiabatic quantum pumping 25 2 6127128] 

• qubit-readout [16129130131) 

• superconducting contacts in equilibrium [§] 

• proximity effect structures [H) 32 33 3 4135) 

• cross-correlations with normal jSS] or superconducting contacts [12137) 

• entangled electron pairs [38139) . 

• phonon counting [40) 

• relation between photon counting and electron counting [41j 

• current biased conductors 02] 

• interaction effects: weak and strong Coulomb blockade [14143144] . 

• multiple Andreev reflections in superconducting contacts [45146) 

Very recently, an important experimental step forward was achieved. Reulet, 
Senzier, and Prober measured for the first time the third cumulant of current 
fluctuations produced by a tunnel junction 07] . Surprisingly the measured volt- 
age dependence deviated from the expected voltage-independent third cumulant 
of a simple tunnel contact [2123) . A subsequent theoretical explanation is that 
the third cumulant is in fact susceptible to environmental effects 02]. This ex- 
periment has already triggered some theoretical activity [24149150) . 

2 Full Counting Statistics 

The fundamental quantity of interest in quantum transport is the probability 
distribution 

P t0 (N 1 ,N 2 ,...,N M ) = P(N), (1) 
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which denotes for a M-terminal conductor the probability that during a cer- 
tain period of time to N\ charges enter through terminal 1, N2 charges enter 
through terminal 2, . . ., and Nm charges enter through terminal M (negative iVj 
correspond to charges leaving the respective terminal). The same information is 
contained in the cumulant generating function (CGF), defined by 



where we introduced the vector of counting fields x — Oti> X2 ; • • • ,Xn)' The 
normalization condition requires J2n = 1 S(x = 0) = 0. 

Charge conservation We are interested in the long-time limit of the charge 
counting statistics, which means that no extra charges remain inside the conduc- 
tor after the counting interval. If we count only the total number of transfered 
charges, we simply have to consider P(N) — J2n ^£ N a ,NP(x)> or , equivalently, 
to put all counting fields equal S(xi = Xj X2 = Xi ■ • • i Xn — x)- Charge conser- 
vation now means that S(xi = X? X2 = X, • • • , Xn — x) — 0- As a consequence 
the CGF depends only on differences between counting fields. This has the direct 
interpretation, that a difference Xa — X/3 is related to a charge transfer between 
terminal a and (3. In general, this means that we need only M — 1 counting fields 
to describe a M-terminal structure. If one of the counting fields, e. g. xm, has 
been eliminated, the charge transfer into terminal M can be restored from the 
CGF, in which all other \a are equal Xa — Xm- In the special case of a two- 
terminal device, the CGF depends only on \ = Xi — X2- We denote this below 
with S(x)- Later we will see that the CGF's are in general periodic functions of 
X, i- e. S(x + 2tt) — S(x)- This ensures that the total charge transfered is an 
integer multiple of the electron charge e, which makes sense, since we are talking 
about electron transport and want to neglect transient effects. 

However, the interesting question, what the charge of an elementary event is, 
can be answered by FCS. Suppose the a CGF has the property S(x + 1-njn) = 
S'(x)- Direct calculation shows that 



where P n (N) is the distribution S n (x) — S(x/n). The probability distribution 
vanishes for all iV which are not multiples of n, thus the elementary charge 
transfer is in units of ne, where e is the electron charge. This has interest- 
ing consequences in the context of superconductivity, in which multiple charge 
transfers can occur |2 1145146] . or for fractional charge transfer |23j . 

Correlations One commonly addressed question is, if two different events (say 
the charges transfered into terminals a and (3) are independent or not. For 
independent events the probability distributions are separable and we find that 
(N^N l ) = (N*)(Nl). In terms of the CGF this means that the CGF is the sum 




(2) 




(3) 
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of two terms: one which depends only on \a and a second one, which depends 
only on x/3- On the contrary, if the CGF can not be written as such a sum, the 
charge transfers in terminals a and f3 are correlated. 



Special distributions (two terminals) If the elementary events are uncor- 
related, the probability distribution is Poissonian. With the average number of 
events is N we have 

Ppotsson = ^e-^ ~ S( X ) = N (e* - l) . (4) 

In the context of electron transport we encounter this distribution mostly for 
tunnel junctions with an almost negligible transmission probability at low tem- 
peratures. Here N = GxVto/e is simply related to the voltage bias and the 
tunnel conductance. 

As second example we consider the binomial (or Bernoulli) distribution. This 
is obtained if an event occurs with a probability T and the number of tries is 
fixed to No: 

Pbinomal = (Jfj T N (1 - T) N °- N 4-» S(x) = No In [1 + T (e* - l)] . (5) 

In some sense this is the most fundamental distribution in quantum transport: 
it gives the statistics of a voltage biased single channel quantum conductor if we 
identify N = eVt Q /h. 



Special distributions (many terminals) For uncorrelated processes the 
CGF takes the simple form 



(6) 



The resulting distribution is just the product of Poisson distributions, taking 
into account total charge conservation. An important example is a multinomial 
distribution for No independent attempts, which can have different outcomes 
with probabilities T a . It has the form 



S( X )=N \n 



T a (e l 



-1 



(7) 



3 Theoretical approach to full counting statistics 

General theory We will follow here the approach to FCS using the Green's 
function technique 8 . Quantum-mechanically we define the cumulant generating 
function by |8lQIU)l^| 

e s W = (T Jf e-^^ W)I(t) } ) (8) 



Full Counting Statistics 

c 2 ^ m) 

t=0 *S ^ \t I 



5* 



Fig. 1. Keldysh time-ordering contour. 

Here Tr- denotes time ordering along the Keldysh-contour Ck, depicted in 
Fig- Ul The time-dependent field x(i) is defined as ±x for t 6 Ci(2)! i.e. x(i) 
changes sign between the upper and the lower branch of Ck- I(t) is the usual 
operator of the current through a certain cross section. Expansion in the counting 
field yields the cumulants. In the second order we find the 2 nd cumulant as 



C 2 (i ) = f" dt f" dt' (5i{t)5i{t')) . (9) 



dt 

lo Jo 

Higher cumulants yield more complicated expressions. 

Current Correlation Functions The cumulants C n {to) are directly related 
to experimentally accessible quantities like current noise or the third cumulant 
of the current fluctuations. Let us demonstrate the relation for the low-frequency 
current noise, defined by 

Si = 2Af f dr/si(r)Si(0)) , (10) 

J — oo 

where 8I(r) = I(t) — (I) and Af = f max ~ j m in is the frequency band width, 
in which the noise is measured. The factor of 2 enters here to conform to the 
review article H. We now transform in © the integration variables from t, t' to 
T = (t+t')/2,r = t—t'. In the limit to = (Af)^ 1 much larger than the correlation 
time of current-fluctuations, the integral over T can be evaluated and we obtain 
from the desired result Si/2. Similar arguments hold for higher cumulants, 
for which the expression corresponding to|^lare less trivial, however. In Ref. [T7| 
it was noted that C3 depends in an quite unusual way on the frequency band 
measured, i.e. it is proportional to 2f max — f m in, which made it possible to prove 
experimentally that the third cumulant is actually measured. 



Keldysh-Green's Functions So far we have formally defined the CGF quan- 
tum mechanically. The relation to standard quantum-field theory methods is 
made in the following way. We introduce the standard Green's function [SI] in 
the presence of a time-dependent Hamiltonian 

H c (t) = H + ± X (t)I, (11) 
2e 

where the time-dependence is only in the 'counting' field xW- The counting field 
couples to the operator / of the current through a cross section, which intersects 
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the conductor entirely. The single-particle operators corresponding to Hq and I 
are denoted by ho and j. 

Using the matrix notation for the Keldysh-Green's functions, we arrive at 
the equation of motion 

G(t,t'; X ) = S(t-t'). (12) 

Here T3 denotes the third Pauli matrix in the Keldysh space and is a result of 
the unusual time-dependence of the counting field. The relation of the Green's 
function (|12[1 to the CGF (JSJ) is obtained from a diagrammatic expansion in x 
(the calculation is formally equivalent to the calculation of the thermodynamic 
potential in an external field, see e. g. [H]). One obtains the relation [S] 

- ^ J(X) , (13) 
e 

where we have restricted us to a static situation, for which G(t, t) is independent 
of time. Note, that the counting current I(x) should not be confused with the 
standard electrical current, which is actually given by I e i = 1(0). Rather, I(x) 
contains (via an expansion in x) oil current- correlators at once. It nevertheless 
resembles a current in the usual sense. E. g., it follows from Eq. (|12[) that the 
counting current is conserved. 



r)v p. 



njG(t, t; x) 



A simplification In a typical mesoscopic transport problem we can access the 
full counting statistics based on the separation into terminals (or reservoirs) and 
a scattering region. Terminals provide boundary conditions to Green's function 
far away from the scattering region. These are usually determined by external 
current or voltage sources and include material properties like superconductivity. 
Let us now take the following parameterization of the current operator in Eq. l|12[l 

j(x) = (VF(x)) Hm ^- (V x - V x a 3 . (14) 

F(x) is chosen such that it changes from to 1 across a cross section C, which 
intersects the terminal, but is of arbitrary shape. Here we have introduced a 
matrix 03 in the current operator, occurring e. g. in the context of superconduc- 
tivity. We assume that the change from to 1 should occur on a length scale A, 
for which we assume Xp -C A hmp,^o (Fermi wave length Xp, impurity mean 
free path li mp , and coherence length £0 = vf/2A). With this assumption we can 
reduce Eq. I|12|l inside the terminal to its quasiclassical version (see Ref. 51 ) 

r y "I 

v F Vj(x,v F ,t,t,x) = \-i-(VF(x.))v F f K , <?(x,v F ,M,x)J ■ (15) 

Here fjc — t 3 & 3 is the matrix of the current operator and g obeys the normal- 
ization condition g 2 = 1. Other terms can be neglected due to the assumptions 
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we have made for A. The counting field can then be eliminated by the gauge-like 
transformation 

S(x, v F , t, t', x ) = e-^w^/^fx, v F , t, f , oyxn~)fK/2 (16) 

We assume now that the terminal is a diffusive metal of negligible resistance. 
Then the Green's functions are constant in space (except in the vicinity of the 
cross section C) and isotropic in momentum space. Applying the diffusive ap- 
proximation |51| in the terminal leads to a transformed terminal Green's function 

G( x ) = e~ lxfK / 2 G(0y xfK / 2 , (17) 

on the right of the cross section C (where F(x) = 1) with respect to the case 
without counting field. Consequently, the counting field is entirely incorporated 
into a modified boundary condition imposed by the terminal onto the mesoscopic 
system. 

Summary of Theoretical Approach This concludes the theoretical approach 
to counting statistics of mesoscopic transport. Let us briefly summarize the 
scheme to follow. The FCS can be obtained by a slight extension of the usual 
Keldysh-Green's function approach, which is widely employed to treat quantum 
transport problems. Making use of the separation of the mesoscopic structure 
into terminals and a scattering region, the formalism boils down to a very pow- 
erful, but nevertheless simple rule: we have to apply the counting rotation (|17fl 
to a terminal, thus providing new boundary conditions (now depending on the 
counting field x) to the scattering problem. We then proceed 'as usual' and 
calculated the current in the terminal, which again depends on x- Finally the 
counting statistics is obtained from Eq. (IB) . 

4 Two- Terminal contacts 

Tunnel contact To illustrate the theoretical method we first calculate the 
counting statistics of a tunnel junction. As usual the system is described by a 
tunnel Hamiltonian H = Hi + H 2 + Ht, where #1(2) describe the left (right) ter- 
minal and Ht describes the tunneling. The current is calculated in second order 
in the tunneling amplitudes and we obtain I(x) = §f J <iETr (tk [Gi(x), G2] ), 
where Gt is the conductance of the tunnel junction and we have included the 
counting field in Gt. The CGF is (using {d/dx)G 1 { X ) = (i/2) [f K ,Gi(x)]) 

S(X) = ^ lo^' 1 ^ = dETr {Gi(x),G 2 } , (18) 

which is the general expression for the FCS of a tunnel junction. We use the 
pseudo-unitarity fj^ = 1 to write 



S(x) = Nu(e lx - 1) + N 21 (e-^ - 1) , 



(19) 
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where N tj = (t G T /l6e 2 ) J dETi [(1 + f K )&i{l - f K ) Gj ] denotes the average 
number of charges tunnel from i to j. The statistics is therefore a bidirec- 
tional Poisson distribution [23- It is easy to see that the cumulants are C n — 
N12 + {-l) n N 2 i. If either N 21 = or N 12 = we obtain the Schottky limit. Fur- 
thermore, in equilibrium N 12 = N 21 and the FCS is (2G T k B Tto/e 2 ){cos{x) - 1), 
which is non-Gaussian, remarkably. 



General CGF for quantum contacts Using the method presented in the 
previous section, we can find the counting statistics for all conductors, which are 
characterized by a set of transmission coefficients {T n }. Nazarov has shown that 
the transport properties of such a contact are described by a matrix current |5.'i| 

_ e 2 2T n [ Gl ,G 2 ] 
Il2 ~~^l + T n ({G ll G 2 }-2)- (20) 

Here G^) denote the matrix Green's functions on the left and the right of the 
contact. We should emphasize that the matrix form of (|2U[1 is crucial to obtain 
the FCS, since it is valid for any matrix structure of the Green's functions. The 
scalar current is obtained from the matrix current by 

h2 = / dETrf K i 12 . (21) 

To find the FCS, we apply the counting rotation i|17|l to terminal 1, i. e. G\ 
becomes x _ dependent. It turns out that the CGF can then be found generally 
from the relations (|13l) . I|2(JI) . and l|21l) . To integrate l|13l) with respect to Xj we 
need the relations i{d/d X ) Gl ( X ) = [fjcAfr)] and tr[G x (x), {Gi(x), G 2 }»] = 0. 
We find 



^({G 1 ( X ),G 2 }-2) 



(22) 



This is a very important result. It shows that the counting statistics of a large 
class of constrictions can be cast in a common form, independent of the contact 
types. Note, that Eq. I|22l) is just the sum over CGF's of all eigenchanncls. 
Thus, we can obtain the CGF's of all constrictions from a known transmission 
eigenvalue density. These are known for a number of generic contacts (see e.g. |54| 
and Tabled, can be determined numerically, or can be taken from experiment. 
Below we will discuss several illustrative examples for a single channel contacts. 



Normal contacts Consider first a single channel with transmission T between 
two normal reservoirs. They are characterized by occupation factors fi( 2 ) — 
[exp((£' — ^i(2))/fcs'7e) + l]^ 1 (T e is the temperature). We obtain the result |2l(i| 
(see Appendix) 



S{x) = ^J dEln [ l + T ^ E ) ( e * X - !) +T 21 {E) (e-** - 1)] . (23) 
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P (T)[G/G Q ] 




Single channel 


5(T-Ti) 


ln(l-Ti(A-l)/2) 


Diffusive connector 


1 1 

2 7V1 - T 


1 2 

— acosh (A) 
4 


Dirty interface 


1 1 


^2(1 + 4) 




Chaotic cavity 


2 1 


4 In (2 + ^2(1 + A)) 


7T vrvi - r 



Table 1. Characteristic functions of some generic conductors. The transmission eigen- 
value densities are normalized to G/Gq, where Gq = 2e 2 /h is the quantum conduc- 
tance. The third column displays the CGF-density, which determines the CGF via 
S(x) = (toG/4eh) J dEtrs({G 1 ( X ),G 2 }/2). 



Here we introduced the probabilities T\z = Tfi(E) (1 — /2(E)) for a tunneling 
event from 1 to 2 and Tzi(E) for the reverse process. We see that the FCS (for 
each energy) is a trinomial of an electron going from left to right, from right to 
left, or no scattering at all. The counting factors e ±lx — 1 thus correspond to 
single charge transfers from 1 to 2 (2 to 1). 

At zero temperature and ^1—^2 — eV > the argument of the energy 
integral is constant in the interval /ii < E < /12 and we obtain the binomial 
form S(x) — 2e *fr ^ In [l + T (e ix — l)] . Note that for reverse bias ^2 > Mi the 
CGF has the same form, but with a counting factor e~ %x — 1. The prefactor 
denotes the number of attempts M — etoV/h to transfer an electron 1 . If the 
transmission probability is unity the FCS is non-zero only for N — M, which 
therefore constitutes the maximal number of electrons occupying an energy strip 
eV that can be sent through one (spin-degenerate) channel in a time interval £o- 
In equilibrium it follows from Eq. (|23(l that the counting statistics is 

2t k B T e i . 2 / /— . x\ t nA\ 
S(x) = 7 asm (VTsin~J . (24) 

The fluctuations are non-Gaussian, except for T = 1, when S(x) = — tofc ^ Te ' x 2 ■ 



SN-contact The FCS of a contact between a superconductor and a normal 
metal also follows from the general expression Eq. (|22(l . Using the Green's func- 
tions given in the Appendix we find the result 21 



2 

1+ J2 A nq (E)(e^-l) 

q=-2 



(25) 



The coefficients A nq (E) are related to a charge transfer of q x e. For example, 
a term exp(2i%) — 1 corresponds just to an Andreev reflection process, in which 

1 The noninteger values of M(tg) occur due to the quasiclassical approximation A 
more careful treatment reveals that M itself is described by a probability distribution. 
For large M the difference is negligible. 
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two charges are transferee! simultaneously. Explicit expressions for the various 
coefficients are given in Refs. |21I56| . The most interesting regime is that of pure 
Andreev reflection: eV, ksT <C A. Here we obtain 

S(x) = ~J dE\n [1 + R A f + f_ (e i2 * - l) + R A (1 - /+)(1 - /_)) {e^ - l)] , 

where Ra = T 2 /(2 — T) 2 is just the Andreev reflection probability and f± = 
f(±E) denotes the occupation with electrons above(below) the chemical poten- 
tial of the superconductor. For low temperatures ksT e <§; eV <C A, the CGF 
becomes 

5(x) = 2e^ ln[l + Ju(ei2x _ l}] (2?) 

The CGF is now 7r-periodic, which according to Sec. |3 reflects that the charge 
transfer of an elementary event is now 2e, a consequence of Andreev reflection. 
Quite remarkably, the statistics is again a simple binomial distribution. In equi- 
librium, we can adapt the result from Eg. 1241 to find 

S ( X ) = - 2to ^ e ' asin 2 (V^sinx) (far* e [-tt/2, tt/2]) . (28) 

The counting statistics is also non-Gaussian, except for Ra = 1. 



Superconducting Contact Now we turn to a slightly more involved problem: 
a contact between two superconductors biased at a finite voltage V. For eV < 2A 
the transport is dominated by multiple Andreev reflections (MAR). The micro- 
scopic analysis of the average current and the shot noise calculations suggest 
that the current at subgap energies proceeds in "giant" shots, with an effective 
charge q ~ e(l + 2Z\/|eV r |). However, the question of size of the charge transfered 
in an elementary event can only be rigorously resolved by the FCS. The answer 
was given by Cuevas and the author based on a microscopic Green's function 
approach. Independently, Johansson, Samuelsson and Ingerman |46| arrived at 
the same conclusion using a different method. 

Now, what would we like to have? In Sec. |3 we have discussed that one 
can speak of multiple charge transfers if the CGF allows an interpretation in 
terms of elementary events, which are described by counting factors e mx — 1, 
where n denotes the charge transfered in the process. How can we ever hope to 
obtain this from the general formula l|22|l? We have to calculate the determinant 
of a 4x4-matrix, which can give only factors of the type e l2x or even smaller 
charges. The answer to this puzzle is that we have to reinterpret the matrix 
structure in (|22J) . since the Green's functions of superconductors at a finite bias 
voltage are essentially nonlocal in energy. The general result for the CGF can 
be written as S(x) = TrlnQ, where Tr= f* a dttr and Q(t) = 1 + (T/4)({Gi f 
G 2 }-2)(t, t). Here Gi®G 2 (t, t') = J dtf'G^t, t"G 2 (t", t'). Let us set the chemical 
potential of the right electrode to zero and represent the Green's functions by 
Gi(M') = e ieVtf3 G s (t - t')e~ ieVt '^ and G 2 (M') = G s (t - t'). Here, we have 
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not included the dc part of the phase, since it can be shown that it drops from 
the expression of the dc FCS at finite bias. The Fourier transform leads to 
a representation of the form G(E,E') = J2 n Go,n(E)6(E — E' + neV), where 
n = 0, ±2. Restricting the fundamental energy interval to E — E 1 € [0, eV] 
we can represent the convolution as matrix product, i.e. (G\ <g> G2)(E,E') — > 
(6iG 2 )„, m ( J B, J B / ) = J2k( G i)™A E i E ')(G2)k,m(E,E'). The trace in this new 

representation is written as J Q e dEJ2 n Trln (Q) In this way, the functional 
convolution is reduced to matrix algebra for the infinite-dimensional matrix Q. 
From these arguments it is clear that the statistics is a multinomial distribution 
of multiple charge transfers: 



S(X) = 



to 



dE In 



1 



OO 

£ 

n— — oo 



Pn(E,V) (t 



1) 



(29) 



General expressions for the probabilities P(E, V) have been derived in Ref. |45j . 

Here we will pursue a different path and study a toy model. Let us neglect 
all set f R - A (\E\ < A) = 1, g R( - A *>(\E\ > A) = ±1, and equal to zero otherwise. 
Physically, this means that we neglect Andreev reflections above the gap and 
replace the quasiparticle density of states by a constant \E\ > A. This simplifies 
the calculation a lot, since the matrix trace now becomes finite. Let us for ex- 
ample consider a voltage eV = 2 Aj 4. In that case, we consider the determinant 
of the matrix 



dct 



2 



(Q-ix) 



V 





1 



1 

e-*> 



\ 



(30) 



where Q±(x) describe quasiparticle emission (injection) and off-diagonal pairs 
e ±x are associated with Andreev reflection. Evaluating the determinant we find 
S( x ) = 4£ In [1 + P 5 (e m * - 1)] , where P 5 = T 5 /{16 - 20T + 5T 2 ) 2 . This ex- 
pression describes binomial transfers of 5 charges with probability P5 . For general 
subharmonic voltages 2A/(n — 1) we find 



s(x) 



2At n 



(n- l)h 

where the probabilities are given by 



In [1 + P n (e m x - 1)] , 



P^ = 



{2-T} 2 ' ^ 3 (4-3T) 2 ' 1 4 (8-8T+T' 2 y^ ' 

p _ T 6 p _ T 7 

r ° ~ (2-T) 2 (16-16T+T 2 ) 2 ' 1 7 (64-112T+56T ;i -7T ;i ) 2 



(16-20T+5T 2 ) 



(31) 



(32) 



Note the limiting cases of these probabilities P n ~ T Tl /4 Tl_1 for T <C 1 and 
P n = 1 for T — > 1. We conclude this section by saying that the general results 
for the CGF [IS] allow for a fast and efficient calculation of all dc-transport 
properties of contacts between superconductors (which may contain magnetic 
impurities, phonon broadening or other imperfections). 
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5 Quantum Noise in Diffusive SN-Structures 

In this section, we illustrate a further advantage of the Keldysh-Green's func- 
tions approach to counting statistics. We consider a normal metallic diffusive 
wire connected on one end to a normal metal reservoir and on the other side 
to a superconductor. The wire is supposed to have a mean free path I ^> Af, a 
corresponding diffusion coefficient D — VfI/3, and a length L. For eV, ksT <C A 
the transport occurs through Andreev reflection at the interface to the supercon- 
ductor. This system shows a quite remarkable property, which is the so-called 
reentrance effect of the conductance. The energy difference 2E of electron-hole 
pairs leads to a dephasing on a length scale £e = \/D/2E. This has the con- 
sequence that the (otherwise) normal wire becomes partially superconducting 
and the conductance increases with decreasing energy. However, once the coher- 
ence length £e reaches L the conductance decreases again. Finally for E = 
the conductance is exactly equal to the conductance in the normal state. This 
is the reentrance effect occurring at an energy of the order of the Thouless en- 
ergy E c = hD/L 2 . In Fig. |21 (left panel, dotted curve) the resulting differential 
conductance at zero temperature is plotted. 

The transport in this system is described by a matrix diffusion equation for 
the Keldysh Green's functions, the so-called Usadel equation 

-— VI = [-iEr 3 ,G] , I = -aGVG. (33) 

In these equations a = 2e 2 N D is the conductivity. The boundary conditions for 
this equation are that the Green's functions in the terminal approach the bulk 
solution Gn or Gs, respectively. This equation is in general difficult to solve, 
even if one is interested in the average current only. However, we can calculate 
the noise and the counting statistics using the recipe outlined in Sec. |3| and 
obtain the noise in the full parameter range of Eq. (|33|l . 

Before considering Eq. I|33l) in its full generality, we consider the limiting 
cases of low and high energies (compared to E c ). For E = the r.h.s. is absent 
and the system is completely analogous to a diffusive connector as discussed in 
0] From Tableland using the eigenvalues (|52J) we find 

S(X) = ^ J d^acosh 2 [2 (f + f-(e 2 ^ - 1) + (1 - / + )(1 - f-)^ - 1)) - l] . 

(34) 

This result shows, once again, that the charges are transfered in pairs. It is 
interesting to compare with the CGF for a diffusive wire between two normal 
metals, for which we obtain 0E1 

Six) = 5j§ / ^acosh 2 [2 (A(l - /aXe* - 1) + f 2 (l - / x )( e -^ - 1)) - l] . 

(35) 

We see that the only difference in the CGF between the SN- and the NN-case 
is in the counting factors, and a prefactor 1/4. Note, that this coincidence only 
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Fig. 2. Noise in diffusive SN-systems. Left panel a): the differential conductance and the 
noise show a reentrant behavior. The effective charge, defined as q e ff(E) — (3/2)dS/dI 
reveals that the correlated Andreev pair transport suppresses the noise below the un- 
corrected Boltzmann-Langevin result 2e. Right panels b) and c): Effective charge of 
the Andreev interferometer shown in the inset (realized experimentally in Ref. |33p . 
The upper panel b) shows the theoretical predictions and the lower panel c) the ex- 
perimental results. The theoretical results contain no fitting parameter (the Thouless 
energy E c = 30/ieV was extracted from the sample geometry and the experimental tem- 
perature of 43mK was included in the calculation). Therefore, it is reasonable that the 
deviations between experimental and theoretical results comes from possible heating 
effects in the experiment, which are not accounted for in the theoretical calculation. 



occurs for the diffusive connector, but is by no means a general rule. At zero 
temperature the results simplify and we find 

S SN (x) = ^ NN (2 X ) , S^(x) = ^acosh 2 (2e«* - l) , (36) 

a surprising simple relation between the CGF for the Andreev wire and the 
normal diffusive wire. It is easy to see that the cumulants obey the general 
relation G„ = 2 n_1 C^ 1N . We observe that we can read off the effective charge 
from the ratio C^ N /C^ N = (^e/z/e)"" 1 and, indeed, find q e ff — 2e. This result 
for the effective charge is a special property of the diffusive connector. 

At energies large compared to E c it is also possible to find the CGF for the 
Andreev wire in general. Then the proximity effect in the wire is absent and 
it turns out that the wire can be effectively mapped on a normal circuit, 
consisting of two identical wires in series to which twice the voltage is applied 
and twice the counting field. Thus, for E 3> E c we obtain S SN (x) from S NN (x) 
by the replacement x ~> 2% and G — ► G/2, which exactly brings us to Eq. 13^1 
and shows that the counting statistics is again the same in the incoherent limit. 

The full quantum-mechanical calculation of the energy-dependent shot noise 
can be performed on the basis of the approach of Sec. 01 lQJ- We expand up 
to linear order in i.e. G{x) — Go — i(x/2)Gi and I(x) = Iq — i(x/2)Ji. 
Substituting in l|33|) we find 

-jrh = [-iEH, Gi] , h = -a (G ^-Gi + Gi— G ) . (37) 
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The boundary conditions at the reservoirs read Gi(0) = [tk,Gl] at the left end 
and Gi(L) — at the right end. Finally the noise is Sj — —e J eLETrfK-Zi(x). 
By taking the trace of Eq. (|37|l multiplied with tk it follows that it does not 
matter, where the noise is evaluated, as it should be. From these equations the 
generalization of the Boltzmann-Langevin equation to superconductors can be 
derived which allows for a faster numerical solution. The results for the 
energy dependent noise is shown in the left panel Fig. |21 A direct comparison of 
the differential shot noise and the differential conductance (for zero temperature) 
shows the difference in the energy dependence. The effective charge defined as 
q e ff = (3/2)dS/dI displays the clear deviation of the quantum noise from the 
Boltzmann-Langevin result of 2e. At energies below the Thouless energy E c the 
effective charge is suppressed below 2e. This shows that the correlated Andreev 
pair transport suppresses the noise below the uncorrelated Boltzmann-Langevin 
result. 

To experimentally probe the pair correlations in diffusive superconductor- 
normal metal-heterostructures it is most convenient to use an Andreev interfer- 
ometer. An example is shown in the left part of Fig. A diffusive wire connected 
to a normal terminal is split into two parts, which are connected to two different 
points of a superconducting terminal. By passing a magnetic flux through the 
loop one can effectively vary the phase difference between the two connections 
to the superconductor. Such a structure has been experimentally realized by the 
Yale group |33|. In Fig. [3 we present a direct comparison between our theoretical 
predictions and the experimentally obtained effective charge. Note, that we have 
included the experimental temperature in the theoretical modeling. The finite 
temperature explains the strong decrease of the effective charge in the regime 
|eV| < fcgT, where the noise is fixed by the fluctuation-dissipation theorem. The 
disagreement between theory and experiment in this regime stems solely from 
differences in the measured temperature-dependent conductance from the theo- 
retical prediction. We attribute this to heating effects. The qualitative agreement 
in the shot-noise regime \eV\ > ksT is satisfactory, if one takes into account, that 
we have no free parameters for the theoretical calculation. Both, experiment and 
theory show a suppression of the effective charge for some finite energy, which is 
of the order of the Thouless energy and depends on flux in a qualitative similar 
manner. Remarkably for half-integer flux the effective charge is completely flat, 
in contrast to what one would expect from circuit arguments based on the con- 
ductance distribution in the fork geometry. Currently we have no explanation 
for this behavior, and therefore more work is needed in this direction. 

6 Multi Terminal Circuits 

In circuits with more than two terminals it is of particular interest to study non- 
local correlations of currents in different terminals. For that purpose we need a 
slight extension of the theoretical approach of Sec. 01 suitable for multi-terminal 
circuits. We will now introduce this method. 
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Circuit Theory To study transport in general mesoscopic multi-terminal struc- 
tures the so-called circuit theory for quantum transport was developed by Nazarov 
53 53] . Its main idea, borrowed from Kirchhoff's classical circuit theory, is to 
represent a mesoscopic device by discrete elements, which resemble the known 
elements of electrical transport. We briefly repeat the essentials of the circuit 
theory. Topologically, one distinguishes three elements: terminals, nodes and 
connectors. Terminals are the connections to the external voltage or current 
sources and provide boundary conditions, specifying externally applied voltages, 
currents or phase differences in the case of superconductors. The actual circuit 
is represented by a network of nodes and connectors, the first determining the 
approximate layout and the second describing the connections between different 
nodes, respectively. 

To describe quantum effects it is necessary to represent the variables de- 
scribing a node by matrix Green's function G, which can be either Nambu 
or Keldysh matrices, or a combination thereof. Consequently, we describe the 
current through a connector by a matrix current I, which relates the fluxes 
of all elements of G on neighboring nodes. The current has been derived by 
Nazarov and is given by Eq. I|20(l for a connector, characterized by a set 
of transmission coefficients {T„}. Note that the electrical current is obtained 
from I\2 = J dETifKli2- The boundary condition are given in terms of fixed 
matrix Green's functions Gi, which are determined by the applied potential, the 
temperature, the type of lead, and a counting field \i- 

Once the network is determined and all connectors are specified, the transport 
properties can be found by means of the following circuit rules. We associate an 
(unknown) Green's function Gj to each node j. The two rules are 

1. G| = 1 for the Green's functions of all internal nodes j. 

2. the total matrix current in a node is conserved: Iij = 0, where the sum 
goes over all nodes or terminals connected to node j and each matrix current 
is given by (J2DJ. 

Finally, the observable currents into the terminals are given by ij = ■ I+j , 
where the sum runs over all nodes connected to the terminal i. To obtain the 
counting statistics, we finally integrate all currents Ii(x) — (d/dxi)S(x) to find 
the CGF S(x)- 

Multi tunnel junction structure A general expression of S(x) can be ob- 
tained for a system of an arbitrary number of terminals connected to one common 
node by tunnel contacts, see Fig. [3] 3 6112) . At the same time it nicely demon- 
strates the application of the circuit theory rules, presented above. Let us denote 
the unknown Green's function of the central node by G c (x)- The matrix current 
from a terminal a (a = 1, ... , K) into the central node is given by the relation 

L(X) = ^[G c (x),G a ( Xa )] , (38) 

where g a — Gq ^ n T n is the conductance of the respective tunnel junction junc- 
tion, for which we have assumed that all T n <C 1 and g a 3> Gq to avoid Coulomb 



16 Wolfgang Belzig 




Fig. 3. Multi tunnel junction structure: a) general setup with K terminals connected 
to a common node, b) beam splitter setup in which terminal 3 is either a normal metal 
or a superconductor. 



blockade. The Green's function of the central node is determined by matrix cur- 
rent conservation, reading Yla=i^ a ~ E«=i 9aG a , G c ]/2 = 0. Employing the 
normalization condition G 2 C = 1, the solution is 

G C ( X ) = T,Li9 a G a ( Xa ) = (3Q) 

yS^=l 9a9P {G a (Xa), Gp (xp) } 

To find the cumulant-generating function (CGF) S(x) we integrate the equations 
dS( X )/d Xa = (-*to/4e 2 ) / dETrf K I a ( X ) EH We obtain 

S(X) = ^J dET^Y^.^^P {G a (x a ),Gp( X p)} ■ (40) 

This is the general result for an M-terminal geometry in which all terminals 
arc tunnel- coupled to a common node. It is valid for arbitrary combinations 
of normal metal and superconductor, fully accounting for the proximity effect. 
Note, that we have dropped the normalization of S( X ) to write the expression 
more compact. 



Normal metals If all terminals are normal metals, the matrices in Eq. I|40|l are 
all diagonal and trace is trivial. We obtain 



S(X) = f- 2 \ dE gl + ]T g a gME)(l - ME)) (jbc*-*,) - l) (41) 

where f a is the occupation function of terminal a. Here, we introduced the 
abbreviation gz = Y]^_-, g a for the sum of all conductances. We note, that the 
statistics is essentially non-Poissonian, despite the fact the we are considering 
tunnel junctions. 

We now restrict us to two terminals (in which case we have to consider only 
one counting field x — Xi ~ Xi)- For zero temperature and voltage bias V the 
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CGF reads then 

S(x) = ^^/gl + ^gmie^-l), (42) 

the result for a double tunnel junction first obtained by de Jong ^3] using 
a master equation approach. We obtain as limiting cases for an asymmetric 
junction (either <?i <C 32 or g\ ^> g%) Poisson statistics Six) — (to ^3132/(31 + 
ffa))(exp(ix) - 1). 

Next we consider a three terminal structure, which is voltage biased such 
that the mean current ^3 in lead 3 vanishes (voltage probe) and a transport 
current I = 31.92/(31 + 32) V flows between terminals 1 and 2. The CGF is [23 

Six) = ^ (92^9% + 43331 (e" ixi - 1) + 40ifl2(e*»-*xi - 1) 

+31/9I; + 4. 9 332(e i x 2 - 1) + 4. 9l52 (e i X2- 4 xi - 1)^ . (43) 

It is interesting to note that the presence of the voltage probe makes the CGF 
asymmetric under the transformation g\ «-> 52, whereas the current is symmetric. 
In certain limits in which the square roots in Ea. 1431 can be expanded one is able 
to find the counting statistics. E. g in the strong-coupling limit g% 3> (31 + 32) 
we find 

Six) = N [e~ lXl + e lX2 - 2] . (44) 

The CGF is simply the sum of two Poisson distribution, demonstrating dras- 
tically the effect of the voltage probe. It completely suppresses the correlation 
between electrons entering and leaving the central node. 

Another interesting geometry is a beam splitter configuration, in which a 
voltage bias is applied between one terminal and the other two. We find 

S N ixu Xa) = ^p\M + 5133 (e* 1 - 1) + 3332 - 1) . (45) 

In the limit that gi + g2 and g$ are very different, we can expand the CGF 
and find for the CGF Six) — Nxe %xi + N 2 e lx , i. e., the tunneling processes into 
the two terminals are uncorrelated. The corresponding probability distribution 
is simply the product of two Poisson distributions. 

SN-contact We now consider the case of a double tunnel junction, in which 
one of the terminals is superconducting. From the general result <|4(J|) and l|52|l 
we find after some algebra 

Six) = ^{al + g\ + \/(3? + glf + ±g\g\{^ i). (46) 

Remarkably, the statistics is fundamentally different from the corresponding 
normal case (|42H . Still, the elementary events are transfers of pairs of electrons, 
which, however are correlated in a more complicated way than normal electrons. 
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If the junction is very asymmetric, the FCS reduces to Poissonian transfer of 
electron pairs. This is similar to the effect of decoherence between electrons and 
holes for energies of the order of the Thouless energy j.'i2j . 

For the beam splitter configuration we are also able to find the FCS analyt- 
ically. The CGF is [E] 

S( Xl , Xa ) = ^x (47) 

^9 2 S +^9s + 4<?I«7i(e l2xi - 1) + 4. 9 f ff 2 (e^ - 1) + 8g 2 gi g 2 (e^+ x ^ -^), 

where we abbreviated g| = g\ + [g\ + g 2 ) 2 • From this result we see that the 
elementary processes are now double charge transfers to either terminal of a 
splitting of a Cooper pair among the two terminal. It is interesting to note, that, 
if we assume that gi + g 2 and 53 are very different (but g± « g 2 ), we obtain 
non-separable statistics 

S(x) = N n e' 2xi + N 22 e l2x2 + N 12 e l(xi+X2) . (48) 

This expression can not be written as a sum of two independent terms. Fur- 
thermore, the last term is positive, which implies that current crosscorrelation 
5*12 = - fa 2 /to){d 2 /dxidx2)S(xi, X2)\xia2^ are positive. Eq. lO provides a 
simple explanation for this surprising effect: it is a consequence of independence 
of the different events, contributing to the current. This result, in fact, holds for 
a large class of superconducting beam splitters [34137160161] . 

7 Conclusion 

We have tried to give a pedagogical introduction to the field of counting statistics. 
Many technical details were left out, but we have tried to cover the essence of 
the derivation and concentrated on looking at concrete examples. For a more 
thorough study we recommend the recent book Quantum Noise in Mesoscopic 
Physics or the original literature. While a number of aspects have already 
been explored, many open questions remain, e. g., experimental strategies to 
measure FCS, strongly interacting systems, or spin-dependent problems. For 
the future, we expect even more activity in the field and, consequently, even 
more interesting results will emerge. 
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8 Appendix 

We summarize here the matrix-Green's function for superconducting and normal 
contact, as they were used in the text. The time-dependent Green's functions are 
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expressed by their Fourier transforms G (t - t') = J(dE/2n) e~ iE ^ t ">G (E). 
The energy-dependent Green's functions in the KeldyshxNambu-space have the 
form 

f (A-R)f + R (A — R)f 
\{A-R){l-f) (R — A)f + A 

where the advanced, retarded and occupation Nambu matrices are 

The phase tp of the superconducting order parameter as well as the electrical 
potential eV enter via the gauge transformation G(t, t') = U(t)Go{t — t')U^(t'). 
Here U{t) = exp [i(f}(t)f 3 /2], where <p{t) = ip + eVt. 

In the calculation of the FCS of contacts between normal metals and su- 
perconductors we frequently need the eigenvalues of anticommutators of two 
Green's functions. For two normal metals {GVi(x); Gn2}/2 is diagonal and the 
eigenvalue is 

[1 + 2h{E) (1 - h{E)) (e* - 1) + 2f 2 (E) (1 - fi(E)) (e"* - l)] , (51) 

for the electron block and the same expression with E — ► —E for the 'hole'-block 
in Nambu space. 

In the case of Andreev reflection, i. e. for eV,ksT e i -C A, we find for 
{Gat(x), Gs}/2 the two eigenvalues 

±VfN(E)f N (-E) (1 - e«x) + (1 - f N (E))(l - f N (-E)) (1 - F^). (52) 
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